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We present an analytic theory of quantum interference and Anderson localization in the quantum 
kicked rotor (QKR). The behavior of the system is known to depend sensitively on the value of 
its effective Planck's constant h. We here show that for rational values of h/iAir) — p/q, it bears 
similarity to a disordered metallic ring of circumference q and threaded by an Aharonov-Bohm flux. 
Building on that correspondence, we obtain quantitative results for the time-dependent behavior of 
the QKR kinetic energy, E{t) (this is an observable which sensitively probes the system's localization 
properties). For values of q smaller than the localization length ^, we obtain scaling E{i) ~ AP, 
where A = 2n/q is the quasi-energy level spacing on the ring. This scaling is indicative of a long time 
dynamics that is neither localized nor diffusive. For larger values q S> Ci the functions E{i) — >■ ^'^ 
saturates (up to exponentially small corrections ~ exp(— g/^)), thus reflecting essentially localized 
behavior. 

PACS numbers: 05.45.Mt, 71.23.-k, 71.55.Jv 

I. INTRODUCTION 

The kicked rotor, or standard map [1] is one of the most prominent model systems of nonlinear dynamics. In 
spite of its very simple construction - the rotor Hamiltonian depends on just a single parameter, the dimensionless 
kicking strength, if - it shows complex physical behavior. Specifically, the classical rotor undergoes a transition 
from integrable to chaotic dynamics as the kicking strength is increased. Quantization in the chaotic regime leads 
to a quantum system that bears strong similarity to a disordered multi-channel quantum wire. Much like in these 
systems, the dynamics of the quantum kicked rotor (QKR) is governed by mechanisms of quantum interference and 
localization [5| |3] . The applied interest in researching these structures has increased since the mid nineties [3] when 
concrete realizations of the rotor in atom optics settings became an option [5] . 

The analogies to condensed matter systems notwithstanding, there are also a number of important differences. The 
QKR is not a genuinely disordered system. Rather, its effective stochasticity roots in mechanisms of incommensura- 
bility. The differences from generic disorder manifest themselves in anomalies showing as Planck's constant assumes 
peculiar "resonant" values [2j|6H9]. (In this paper, we will not discuss the role of classical anomalies, such as accelera- 
tor mode formation [TO] at special values of K.) At these values, the system becomes effectively periodic (in the space 
of angular momentum quantum numbers) and generally escapes localization. Various aspects of QKR resonances 
have been discussed before in terms of general observations, phenomenological reasoning and direct diagonalization 
(see Ref. [8] for review.) However, a microscopic theory of the interplay of resonances and localization has not yet 
been formulated. The construction of such a theory is the subject of the present paper. 

Below, we will introduce an analytic and parametrically controlled theory of localization phenomena in the QKR 
both on and off resonant values of Planck's constant, h. It turns out that at resonant values, h = Anp/q, with p, q 
coprime numbers, localization phenomena in the QKR can be described in terms of a field theory defined on a ring 
subject to an Aharonov-Bohm flux (the latter implementing the Bloch phases pervasive in the physics of periodically 
extended systems.) We will apply perturbative and non-perturbative methods to analyze this field theory in terms 
of different observables. Emphasis will be put on the discussion of the time dependent increase in the rotor's kinetic 
energy, an observable that carries immediate information on the transport characteristics of the system. 

Before turning to a qualitative discussion of results, it is worthwhile to put the present theory into the context of 
earlier work. The connections to disordered systems were first discussed in Ref. [TT] . That work introduced a mapping 
of the QKR Hamiltonian onto an effectively disordered tight binding system. However, the long-ranged correlated 
'disorder' inherent to that model made it difficult to treat analytically. A different approach was introduced in Ref. |12j . 
where a diagrammatic approach conceptually similar to the impurity diagram language of disordered systems was 
introduced. In this work, perturbative quantum interference processes ("weak localization") were discussed in terms 
of low energy fluctuations in angular momentum space (for a refined implementation of the diagrammatic approach, 
see Ref. fTT.) These results inspired a subsequent field theoretical study fT^ in which the large scale dynamics of the 
quantum kicked rotor was mapped onto a one-dimensional supersymmetric nonlinear a-model (otherwise known as 
the effective field theory of Anderson localization in disordered quasi-one-dimensional wires |15).) 

None of these works payed attention to the resonance phenomena at rational /i/(47r), a lack of resolution that has 
prompted some criticism [51 116] . Specifically, it has been argued that the insensitivity of earlier work to the algebraic 
characteristics of h might be indicative of ad hoc assumptions hidden in the construction. One of the motivations 
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behind the present work is to show that these objections have been ill founded. Rather, we will see that the field 
theory below - by and large a refined variant of the earlier theory [14] - is capable of describing the system from large 
scales down to the fine structure of individual quasi-energy levels. 

Nonetheless, it is worthwhile to point out that our theory is far from 'mathematically rigorous'. Rather, it is 
based on a mapping of QKR correlation functions onto a functional integral (exact) followed by an identification 
of field configurations of least action, plus an integration over fluctuations. The last two steps are approximate, 
if paranietrically controlled (by a small parameter h/K.) In this sense, our theory below falls into the general 
framework of semiclassical approaches to nonlinear systems. The semiclassical limit enables us to explicitly describe 
system observables which more rigorous lines of reasoning |17j can only characterize in implicit terms. We also 
emphasize that no averaging procedures other than a global average over the quasi-energy spectrum will be performed 
throughout. 

The rest of this paper is organized as follows: we start with a qualitative discussion of the system and a summary 
of our main results (Section |ll]). In Section III we introduce the functional integral approach to the problem, which 
will be reduced to an effective semiclassical act ion in Section |IV[ The low energy physic s of t his action will then be 
discussed in Sections [v| (off-resonance) and VI (on resonance). We summarize in Section VII Some technical details 
are deferred to Appendices |^][C] 



II. QUALITATIVE DISCUSSIONS 

The QKR is a periodically driven ( "kicked" ) quantum particle moving on a circle. In units where the particle mass 
and the circular radius are set to unity, its Hamiltonian reads 

P 

H{t) = - + kcose^S{t~^mT), (1) 

m 

where T and k are the period and the kicking strength, respectively. The angular operator, 0, and the angular 
momentum operator, I obey the canonical quantization condition [0,1] = ih where I = ifi-^ ■ In rescaled dimensionless 
time, t — > t/T = t, the dynamics of the system is determined by two parameters: the effective Planck's constant h — hT 
and the classical nonlinear parameter K = kT . (In the chosen units, both parameters are dimensionless.) 



A. Symmetries 

At integer times f = 0, 1, 2, . . . , the solution of the Schrodinger equation ihdi\4'{i)) = T'^H\t/j{i)) can be expressed 
as IV'(^)) = C^*lV'(0))i where U is the Floquet operator, i.e. the unitary operator governing the time evolution during 
one elementary time step: 

(ihn'^ \ I iK cos §\ ( ififi^ \ 
j exp j exp j , (2) 

where we have introduced n = l/h as a rescaled angular momentum operator with integer-valued spectrum quantum. 

For Planck's constants commensurable with 47r, h = Airp/q, with coprime p, q € N, the Floquet operator is invariant 
under the shift n ^ fi + q, or 

[f7,f,]=0, (3) 

where Tq is the translation operator by q steps, Tq\n) = \n + q) on eigenstates h\n) = n\n). This discrete translational 
symmetry implies that the (quasi)energy states of the Floquet operator can be expanded in a basis of Bloch states, 

where 9 € [0, 27r/g] and ipn+q,d = V'n.e- The infinite extension of Bloch states in angular momentum space means the 
absence of localization at these "resonant values" of h. ^ [BH5] (However, for periodicity intervals q much larger than 
the intrinsic localization length of the system, the wave function amplitudes at the boundaries of the "unit cell" are 
exponentially small, which means that the system behaves effectively localized.) Below, we will explore the interplay 
of localization and Bloch periodicity in some detail. Our discussion will include the case of irrational h/Ai: as a limit 
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of co-prime values p/q with diverging q. We will assume that q ^ K/h = £, which means that the angular momentum 
unit cell is larger than the length scales at which diffusive dynamics begins to form. 

In dirty metals, the manifestations of localization depend sensitively on the presence or absence of time reversal, 
t — >■ —t,r r,p ^ —p (where r and p are coordinate and momentum respectively.) In the context of the kicked 
rotor, momentum and coordinate space change their roles, which is why time reversal followed by space inversion, 
Tc : t — > —t, 9 — > —9, 1 — > I becomes a relevant symmetry. It has been shown [121 1131 HI] that Tc plays a role 
analogous to time reversal in metals. The Hamiltonian ([T]) is Tc-invariant, and in our analysis below, we need to 
take this symmetry into account. (Tc~invariance may be broken e.g. by a shift cos(0) — )■ cos{9 + a), a — const. 
However, in this paper, we focus on the invariant system.) In a basis of angular momentum eigenstates, {1?^)}, 
l\n) — hn\n), n S Z , the symmetry Tc acts by matrix transposition Ann' ^ {A-'^)nn' = ^n'n (much like conventional 
time reversal acts by transposition in a real space basis.) 



B. Off resonance 



Let us temporarily assume that h e 47rIR\Q is irrational (or rational with q ^ ^, so that the periodicity of the 
system plays no role.) Throughout this paper, we consider the system prepared in a pure initial state with zero angular 
momentum, i.e., |0)(0|. Under the Floquet dynamics it will evolve into (?7)*|0) (0|(i/^)'. Its (squared) deviation from 
the point of departure is measured by the expectation value E{i) = itr(n^(i/)*|0)(0|(C/^)*). Comparison with the 
Hamiltonian ([T]) shows that this quantity may be interpreted in terms of the system's change in kinetic energy. In 
the classical limit, E{t) value grows indefinitely in time, a phenomenon that may be interpreted as a mechanism of 
'heating'. More precisely, 

E{i) ^ Dot, 

where Dq = {K/2h)'^ has the status of a diffusion constant in angular momentum space. Notice that the energy 
expectation value can be represented as E{i) = ^ n^K{n, 0; i), where K{n, n'; t) is the density-density correlation 
function in angular momentum space (for a precise definition, see Section HI below.) The linear growth in time is 
indicative of diffusive behavior K{q,ui) ^ {—iuj + D^q^)^^, where {uj,q) is Fourier conjugate to (t,n). Empirically, it 
has long been known [3j that in the quantum model, the linear growth in energy comes to an end at times t^ ^ Dq. 
Simple scaling arguments show that this saturation implies a crossover Do — !■ D{uj) ^ iu) at frequencies to ^ -Dq"^. 
The scaling D{uj) ^ iuj is & manifestation of quantum localization ^15, . 

A microscopic quantum theory of the system must be capable of establishing diffusive dynamics at intermediate 
time scales, and to describe the quantum interference processes rendering the diffusion constant frequency dependent. 



C. On resonance 



We now discuss what happens for rational values h = Airp/q. The commutativity ([s]) means that we are considering 
a variant of a g-periodic quantum system. For the convenience of the reader, a few basic structures of periodically 
extended quantum systems are recapitulated in Fig. [l] The Bloch function basis can be understood as a basis whose 
elements exhibit a definite phase change 6q e [0, 27r] across a single unit cell (top panel.) For a given 9, one may 
then consider the system within the reduced scheme of just a single unit cell, subject to twisted boundary conditions 
'(/'(O) — il^{q) eyi\)[i9q) (middle left). Equivalently, one may interpret the system in terms of an Aharonov-Bohm ring 
(middle right) subject to a gauge flux 9q. The q energy levels ej{9) of this system (bottom) form a.2n/q periodic 
family, each j defines a Bloch band. 

Pioneering work on the physics of the kicked rotor at these "resonant" values of Planck's constant has been done by 



Izrailev and Shepelyansky (IS) |S] (see Ref. [S] and Section VIA below for review.) In this work, it has been argued 
that at resonance, the energy stored in the system increases as 

E{i) = T]P + ai+b{i), (4) 

(designation of constants rj, a and function b{t) taken from the review Ref. [S]). Let us briefly discuss the meaning of 
this expansion, as discussed in IS. Eq. Q has the status of a large time asymptotic. For times large enough such that 
is smaller than the level spacing of the unit cell, A = 27r/g, we are in a "deep quantum regime", where transitions 
between individual Bloch bands can be neglected. In this regime, the energy increase is governed by expressions such 
as (symbolic notation) 
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FIG. 1: Wave functions and spectrum of the periodically extended system. 



Thus, the contribution is a measure of the Bloch-phase auto-correlation of individual levels. The auto-correlation 
of levels also yields a quantum correction linear in t, and a residual contribution 6(t) of indefinite time dependence. 
The coefficients of both terms, ry and a, reflect the dispersion of Bloch bands, and vanish if bands are flat. Finally, 
in regimes with localization, ^ <C g we expect the increase in n? to come to an end at large times i ~ ^^/Dq. At 
larger times, E{t) will saturate. (With exponentially small corrections ~ exp(— g/^) proportional to the wave function 
overlap across the ring.) 

The construction of IS that led to Q has been formal in that it relied on explicit knowledge of eigenfunctions and 
-values of the Floquet operator; apart from a few exceptional values of g, no concrete results for the coefficients of the 
expansion have been available. Equally important, at times shorter than mechanisms of energy diffusion different 
from those discussed in IS play a role. Specifically, for small values of time, the finite extension of the unit cell is not 
yet felt, and angular momentum will diffuse as in an infinite system. This generates a (leading) ~ t dependence in 
h{t), which is purely classical. (Within the framework of a quantum approach, diffusion has to do with inter-band 
transitions.) At short times, the profile of E{i) will, thus, be different from Q. 

Below, we will analytically derive the function E{t) for periodicity volumes obeying £ < q < S,. As a result, we 
obtain a universal scaling function 



E{i)^D,qF 



^3 



3' 



where E^ = Dg/q^ is the inverse of th e classical diffusion time through the periodicity volume. Referring for a more 
substantial discussion to Section |VIB| we here merely note that Eq. ([5| predicts diffusive scaling E ~ Dqi for times 



t < A ^. For larger times, we obtain a quadratic increase, E ~ Tyt^, governed by the universal coefficient 77 = Dq/q 



For the discussion of the meaning of this coefficient and of the correction terms in ([5]), we refer to section VI B 



In Section |VI C[ we will briefly discuss the regime of large resonance volumes q > Our results are in line with 
the general expectations formulated above. 
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D. Close to fundamental resonance 

The sensitivity of the rotor to the algebraic properties of h manifests itself not only on resonances, but also in the 
close vicinity of low-order resonant values. To be specific, consider values h — Anp + e, close to the fundamental 
resonance q = I. Straightforward substitution into ([2 1 shows that the Floquet operator defined for the pair {h, K), is 
identical to one with {h^, K^), where [TOj 

he^e, K,^-^. (6) 

This duality phenomenon has been dubbed "e-classics" : a rotor with large value of Planck's constant becomes equiv- 
alent to one with small h ^ e. Faithfulness to this symmetry is a benchmark criterion which the present approach 
obeys. In practice, this means that system characteristics such as the localization length will depend on h and K 
through certain invariant combinations (that have been conjectured phenomenologically before [10].) 
We now turn to the quantitative formulation of our theory. 

III. FUNCTIONAL INTEGRAL FORMULATION 

As in condensed matter physics, localization properties of the QKR can be conveniently probed in terms of "two 
particle Green functions" . In this paper, emphasis will be put on the density correlation function 

+ 00 

KUni,n2)= ({n,\{Ue'^+Y\n2}{n2\{Ue'^-)-''\n,}) , (7) 

t,t'=0 

where the frequency arguments, uj± = ujq±'^, with lu understood as lu + iS, S \i 0, and we have introduced an average 
over the quasienergy spectrum, {■)uo '■— /g"" ■ Summation over obtains 

Kjn,,n2) = ({ni\G+iuj+)\n2){n2\G-{LJ^)\n,)) , (8) 

where 



(w±) = y fe'"±J7) = . , (9) 



and the + (— ) superscript designates the retarded (advanced) Green function. Physically, the function K^{rii,n2) 
describes the probability of propagation n2 — >■ rii in time ^ lo^^ . Exponential decay of {—iujK^) at a; — ?> signals the 
onset of localization. The density correlation function also contains information on the time dependent expectation 



value of the energy. Comparing with the definition given in the beginning of Section II B it is straightforward to show 



^E"' / f^e— *X^(n,0). (10) 

In the present formalism, the consequences of time reversal symmetry are best exposed by "taking the square root" 
of the Floquet operator. 



where 



4 y \ 2h 

We next introduce the matrix Green function 



-I / ihn'^ \ I iK cos 6\ , . 

= exp ±— exp ±— ^ . (11) 



G± = „ " ' . (12) 
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We will refer to the newly introduced two- component structure as "time reversal space", or "T" -space and label it 
by indices t,t' , . . . . The matrix G is related to the original Green function as G = Gn, i.e. 

K^{ni,n2) = {{n^\GtM2){n2\G^i\ni)) . (13) 

As a first step towards the construction of a field integral formulation, we introduce a superfield 



where component structure refers to the matrix structure in (12), and the superfields = {V't^,QAn} carry complex 
commuting (anticommuting) components V'q=i (V'a=2)- The index A will be used to discriminate between the retarded 
(A = +) and the advanced (A — —) sector of the theory. Finally, the commuting variables '!/'ct=i = "00=1 ^^re the complex 
conjugate of 'ipa=i while in the anticommuting sector V'a=2 and ^a=2 are independent variables. 
We next define the compound Green function 

G = E]^^ ® Ibf ® G+(w+) + Ef^ ® Ibf ® G-(w_) , (14) 

where matrices with subscript "AR" ("BF") act in the two-dimensional spaces of A (a) indices. Matrices E^^,X ~ 
AR, BF contain zeroes everywhere except for a unity at position (i,j) and Ix are unit matrices. The correlation 
function K^^ can then be represented as a Gaussian integral 

KUni,n2) - ( J D{^,^)e^p {-^G-'iP) T[^^,^Jj^]) (15) 

where the pre-exponential factor 



The absence of a normalization factor in ( 15 ) follows from the fact that the integration over anticommuting variables 
normalizes the integral to unity. 



( / i?(V;,V')exp(-V;G-V))^ =1- 



Wc next reformulate the field integral in a way that will later enable us to make the consequences of the symmetries 
of the problem manifest. Beginning with time reversal, Tc, we use the symmetry of the Green function, G = G'^ to 
write 



^,^±(^,i;^al^), (16) 



where we defined 

and the Pauli matrix a^p accounts for the anti-commutativity of the variables ipa=2, i-e. V't,2AnV't,2A'ri' = 
— V't,2A'n'V't,2An, and the same for -0^. The newly defined integration variables exhibit the "reality condition" (the 
representative of Tc symmetry in the present construction) 

* = ^-^r, (17) 

where the matrix r is defined through 

r = E]^^ ®a\+ E^p ® (icr|) , (18) 



and matrices with subscript "T" act in the two-component "time reversal" space introduced in (16 1. Armed with 
these definitions, the partition function may now be rewritten as 



-**-2*t 



-2*1 _e 
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The next step in the construction of the field theory is the average over the phase ujq- We will perform this 
average with the help of the "color-flavor" transformation 12 Ij , an exact integral transform stating that for arbitrary 
supervectors, ^fi, \l/2, V&i', ^'2', 



'1^0 J 



where integration domain and algebraic structure of the supermatrices Z and Z will be discussed in a moment. 
Applied to the present context, the transformation obtains 

where the supermatrices Z = {Zatn,a't'n'} are subject to the (convergence generating) condition Zaa = {^)"~^ ^ta- 
The integration measure is given by 

dfi(Z, Z) = d(Z, Z)sdet(l - ZZ) 

= d{Z, Z)cxpstrln(l - ZZ), 

where "sdet" and "str" are superdeterminant and supertrace, respectively Notice that these operations include 
the angular momentum indices of the theory. Finally, the integration domain in the boson-boson sector is restricted 
by the constraint \Zx\Z\t\ < 1 (indices in the BF-sector.) 

Using Eq. (17), the Gaussian integration over variables 5* may be performed to obtain 

Both the integration measure, and the action of this functional integral are invariant under transformations 

Z ^ e i Ze ^ , Z ^ e i Ze ^ . 

We may use this freedom to transform the Floquet operator by e <i to the form 

(ihf?\ fiK cos 0\ 
t/^exp)^— jexp)^^— j, 

which will turn out to be convenient in the following. To keep the notation simple, we denote the unitarily transformed 
operator again by U. 

So far, all operations have been exact. In the rest of the paper, we aim to reduce the functional integral to a 
more manageable form describing long-ranged correlations in the system. To this end, we need to identify field 
configurations (Z, Z) of low action. Field configurations of this type necessarily have to obey the constraint 

Z = T'^Z'^T. (21) 

The reason is that any violation of this symmetry leads to huge mismatch between the three terms in the exponent 
of Eq. ( 20 1 and thereby to large values of the action. 

Let us briefly discuss the conceptual status of the symmetry relation (21 ). The symmetry essentially involves the T- 
degrees of freedom of the theory. Retracing their origin, we notice that the T-indices have been introduced to account 
for the symmetry of the operator ([2]) under transposition - the physical Tc-symmetry. The ensuing doubling of indices 
of Z fields, and the corresponding symmetry (21 1 lead to the appearance of two distinct low-energy field configurations 
in the theory. Formally, these are the field components diagonal and off-diagonal in T-space, respectively. Physically, 
these fields have a status analogous to the diffusion and Cooperon modes in disordered metals. 

The partition function reduced to field configurations obeying (21 1 is given by 

dZe'^^^\ (22) 
S[Z] = -hti- ln(l - ZZ) + ^str ln(l - e^'^U^ZUZ). 
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We next turn to the discussion of the second fundamental symmetry of the problem, the symmetry under q- 
translation in angular momentum space, Tq. The presence of this symmetry suggests to use a basis of Bloch functions 
as a preferential basis. A wave function in angular momentum space would thus be expanded as 



Jb 

oo 

i>n,e^ e-'«to+"Vn, (23) 

j = -oo 

where Jf,{dO) = ^ Jq^^'' dO and the states ipn,e are g-periodic, '4'n+q,e — i'nfi- This representation defines an analogous 
expansion of the fields Z: 



Jb 



jj' 

where a = {a,t) is a container index comprising BF- and T-index, and the fields Zn^e,a-n',e',a' are g-periodic, 

^n+q,9 ,a\n' ^6' ,a' — ,a\n' -\-q,9' ,a' — ^n,Q ^a\n' .9' .a' • 

It may be convenient to think of the fields Zn^9-n',9' (we suppress indices in the notation whenever possible) as 
matrices in a tensor space ^Hb, where K'^ is g-dimensional reduced angular momentum space (the 'unit cell'), and 
■Hb is spanned by Bloch angular wave functions ip{0), also with periodic boundary conditions tp{^) = V'(^ + 27r/g). 
Substituting the above expansion into the action ( 22 ) and using Eq. (§, it is straightforward to verify that the action 
assumes the form 

S[Z] ^ -istr ln(l - ZZ) + ^str ln(l - e^'^ulzU^Z), (25) 



where the supertrace now extends over the Bloch quantum numbers, 

str(X) ^Y ('^^^ 

Jb 



and the residual "str" traces over the internal indices of X. Further, = {U{6)n,n'S{0 — 6')} is an effective Floquet 
operator, diagonal in ?^Bi with matrix elements 

U„^,{e) = e^''#i|]'e^f ™^(^^-+»)+»("'-")^. (26) 

fe=0 

It is instructive to compare this with the matrix elements of the original Floquet operator (defined in unrestricted 
n-space): 

Jo 27r 

The difference is that the compactification of the theory to a unit cell {0, . . . , g} in n-space implies a discretization of 
the angular integration, / |f — ^ ^ X^fc- Second, the quantized "momenta" 27rfc/g get shifted by the Bloch momentum 
9 e [0, 27r/(7]. One may think of U{9)nn' as an effective Floquet operator on a ring (in n-space) of circumference q 
which is threaded by an Aharonov-Bohm fiux 9. (Of course, this fiux is purely fictitious and does not break time 
reversal invariance.) It is straightforward to verify the unitarity of the effective Floquet operator, 

n" 

where angular momentum sums X„ = X]n=i are now all meant to run over the unit cell. Finally, notice that the 
action (25) for the field basis introduced in (24) defines a faithful representation of the theory (22); there are no 
approximations involved. 
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Let us close this section with a brief discussion of the general structure of the theory. To this end, we define the 
matrix 



Q = 9'^ar9 S 



1 Z 
Z 1 







I tan 1 1 ^ 


> + l 







With the formal representation U — exp{iG), the zero frequency action can be represented as (cf. Appendix [A 

S[Q]\^=o = ^strln 

This action is manifestly invariant under transformations Q - 
[go, [/] = 0). The go's fulfilling this condition are given by 

1 Zo 
Zo 1 



(27) 



(28) 



50 = 



3o ^Qdo provided that [go,G] = (or, equivalently. 



(29) 



where Zq = {ZQ^at,a't'6nn'} are unit matrices in angular momentum space. The zero- mode configurations (27 1 define 
the field space of the nonlinear cr-model (of orthogonal symmetry) i 15] . The global invariance under the "zero-mode" 
fluctuations go spanning this manifold bears important consequences for the formulation of the fluctuation expansion. 
It implies that we may organize our later analysis of general fluctuations (fluctuations non-commutative with U) 
around any zero mode reference configuration Qq. Below, we will expand around go = 0, or Q = ctaR' corresponding 
to Zq — 0. In the end of the calculation, the extension to generic Q-field configurations may then be obtained by 
generalization of the Gaussian fluctuation action to an action displaying the full go rotation invariance. 

For later reference, we finally note (cf. Appendix |b| that the Q-matrix representation of the density correlation 
function i^cj("-i, "-2) is given by 

KUni,n2) - /(d0+)(d^?_)e^("--"i)(«+-«-) (str(Q„,,^,„,e_ Ef^<E>V)stY{Qn,e_,n,e^ ^ar® ^)) (30) 



£^BP, the average is over the functional with action (25) and the representation (27) is implied 



where V = E)^^ 

Note that the fields Q, by definition, are also g-periodic, i.e., Qne 
larger than q. 



in+q9 .7l+q 



, and in Eq. (30 1 \ni — may be 



IV. EFFECTIVE ACTION 



In this section we reduce the functional integral formulation (22 1 (or, equivalently, (25)) to a more manageable 
field theory of localization in angular momentum space. In doing so, we need to pay attention to fluctuations Z^n' 
inhomogeneous in angular momentum space. We will begin by looking at the zero frequency action, = 0, which 
describes the physics of infinitely long-lived correlations. 

Conceptually, the field Z^n' describes the dynamical evolution of states \n) (E) {n'\. (For example, the action 
(22) essentially measures the overlap of the one-step evolved state, UZW ^ tj\n) ® {n'\lJ^ with the un-evolved 
configuration.) Configurations of low dynamical action are to be expected if (a) n is close to n' (in which case, 
\n) (8) {n'\ evolves close to a "classical trajectory" through angular momentum space, and (b) variations in n are 
shallow, such that the classical action of the reference trajectory is low. Indeed, in the limiting case Znn' — -B„5„„' 
and Bn independent of n, the operators U and tj^ in (22) cancel out and the action vanishes. Our goal is to compute 



the action associated to soft fluctuations around this zero mode limit. 



Fluctuation action 



As discussed in the end of the foregoing section, it will be sufficient to analyze angular momentum space fluctuations 
around the reference point Z = 0. The structure of the full action then follows from its invariance under uniform 



rotations g^. We thus start by considering the quadratic expansion of the action (251 around Z = and at = 0. 
Denoting the ensuing quadratic action by Sa, we have 



Sfi[Z] = ^str{ZZ - ZU^ZU) 



(31) 



E 



nn' n" n" 



{d0){de')stT{Z„e^„,,g,Z^,g,^,,,,,g){S„,,,,,Sn',n'' ~ U {6%,, ^,,,U {6)1,, J. 
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This is an exact representation of the quadratic fluctuation action. We now need to identify those configurations 
Z whose action is asymptotically vanishing. All other, "massive" field configurations can then be integrated out by 
perturbation theory. Following our discussion in the end of the previous section, we split the fields as 

Z = B^C 

into "massless" configurations B and "massive" fluctuations C . The massless sector will be parameterized as 

B„e,„'e' =(r^Y. (0) <5„„. , (32) 

where is a shorthand for a summation over "mode indices" </> = 0, '2-''^ /q, 47r/(j, . . . , 0o smaller than a certain cutoff 
index 0o to be discussed in a moment. The fields B describe fluctuations diagonal in angular momentum space, 
n = n' . To understand the meaning of this limitation, recall the interpretation of Z„j^„2 as the representative of 
a bilinear |ni)(n2|. Turning to a Wigner representation (symbolic notation), \ni){n2\ — )■ W{n,9) = X^An ^^^^^l*^ + 
An/2)(n — An/2|, we notice that the degree of off-diagonality An = ni — 712 encodes information about the direction of 
propagation in angular momentum space (as measured by the angular variable 6 of the rotating particle.) On physical 
grounds, we expect that sense of direction to decay rapidly if the kicking strength K is sufficiently large (for fixed h). 
Indeed we will see in Appendix [c] that the effect of off-diagonal fluctuations Cn,n' amounts to largely innocent ^ 
corrections to the theory. The essential information is carried by the modes i?„^„ describing fluctuations uniform 
in angular space. The Fourier mode decomposition in terms of an index (/> < 0o states that we will concentrate on 
fluctuations at scales n ^ 27r/0o- Deferring the self-consistent determination of 0o to Appendix [c] we here state that 
(j>o ^ h/ K = £~-^ , i.e. fields of low action will fluctuate on scales ^ £, where the scale £ plays a role analogous to that 
of the elastic mean free path in disordered systems. It is important that the cutoff index cjjQ is invariant under the 
transformation (Note also that the condition 27T/q < 0o ^ h/K implies £, as we assume above. Otherwise, 
the massless sector will contain only the constant mode B^ n = const.) 

For the moment we ignore the contribution of the massive modes and substitute Z B, with B parameterized as 



in (32) into the action. This leads to 

Sfi[B] f{d0){de')stY{Be,g>{^)Be'.A^'))Ke,e'{<l>,<l>'), 
where the integral kernel is given by 



Substituting ( 26 ) , it is a straightforward calculation to reduce this expression to 

q-l 



k=0 

q[^\ {9-9' + <pf5^,.^, +0{Kc^/hf. 



2h 

In the last expression, we have approximated the discrete /c-sum by an integral. (For the small values of (t>,9,9' 
under consideration, this is an innocent operation.) Substituting this result into the action and switching back to a 
momentum space representation B{(j)) — !■ Bin), we obtain 



Sn[B]^^ J\nstT{{~idr, + [9, ])B{+td^-[9, ])B), 



where we have passed from sum to integration because B{n) is interpreted as a smooth function of the variable n. 
Here, we have introduced the (bare) diffusion constant 



2 



i'.Hll, (33) 



and the commutators act as 



i[9,B])e.e> = [9 - 9')Bg,e' 
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The action above was obtained by quadratic expansion around the reference configuration Q = (t\^ (cf. discussion 
in the end of Section III ) To obtain its rotationally invariant generalization, we rewrite the action as 



5fl[Q] 



A) 
' 16 



dnsir{{^idn + [e, ])Q{idn - [9, ])Q), 



(34) 



where Q = d'^ARd ^ ^^^^ 9 ~ ( _i ^ ) ^^"^ 9 ^ — y i j describe infinitesimally small deviations off cr^pj . The 

generalization to an arbitrary configuration Q is then given by 5 — > g and Q ^ Q = g<JAR9^^' where g need no longer 
be close to unity. Wc thus conclude that action ( 34 1 defined for arbitrary 



Q 



9<^ar9 ^ 



1 B 
B 1 



(35) 



is the unique rotationally invariant generalization of the quadratic action above. (Starting from (25), the same action 
can be obtained by the straightforward if more tedious fluctuation expansion around arbitrary Z .) 

To conclude the derivation of Eq. (34), we need to explore its stability with respect to the "massive mode" 



fluctuations C that were so far neglected. This analysis is conceptually straightforward yet tedious in practice. 
Referring to Appendix [C] for the actual formulation of the massive mode integration, we here summarize the main 
conclusions: 



• The soft mode action is stable at length scales An > I = 

• Massive mode fluctuations are damped out at scales < 



K 

h ■ 



• Their principal feedback into the soft mode action is a weak renormalization of the diffusion constant, Z?o 
where the so-called quantum diffusion constant |8l I20j is given by 

~ Do{l ~ 2J2 [Kg] - 2Jl [Kq] + 2Jl {Kq)) + ... 

with 



D. 



and Jn{x) are the Bessel functions of the first kind (cf. Eq. (CI61.) 




(36) 



(37) 



• At certain irrational values [23], h assumes the form of an infinite fraction. At these values, long lived angular 
correlations form, and we are not able to get the massive mode fluctuations under control. At these values, the 
theory considered in this paper does not apply. 



B. Frequency action 

To complete the derivation of the effective action, we need to take the flniteness of the control parameter u into 
account. This will generate a second contribution to the action, Si^[Q], adding to the fluctuation action above. In 
deriving S'(^[(5], we will resort to two (parametrically controlled) approximations: first, we will ignore the massive 
modes in the computation of the finite frequency action [531 Second, we will ignore operators of O{ojd'^). At large 
spatial scales, 9„ < h/K and low frequencies, a; ^ 1, the neglected terms are irrelevant. Technically, the assumptions 



above imply that in the expansion of the action (25 1 to lowest order in w, we will set Z ~ B, and ignore terms arising 
from the non-commutativity of Z and U (which would generate contributions of 0{ujd^).) 
We thus start out from 

Su:[B] ~ -^str ln(l - BB) + ^str ln(l - e"^BB) 




Comparison with the definition ( 35 ) shows that this can be rewritten as 



Combining this with the fluctuation action (34), we arrive at the main result of this paper, the invariant soft mode 
action. 
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S[Q] ^YeJ' {^gi'^n - [i ])Qiidu - [0, ])Q - 2iujQal^^ . (38) 

Notice that the action above is universal in that it depends on system parameters only through the diffusion constant 
Dq = Dq{K, h). The diffusion constant, Dq, and therefore the theory as a whole, are invariant under the transformation 
{h,K) — > ih^^Kf). (Note that, for < 1, the constructions discussed in this paper do not work.) 

Conceptually, S[Q\ describes the physics of a disordered metallic ring (the fields Q obey periodic boundary condi- 
tions Q{n) = Q{n + q)) subject to a fictitious magnetic flux 9. (Indeed, the action is mathematically similar to the 
nonlinear cr-model action of persistent current system, a disordered metallic ring threaded by an Aharonov-Bohm 
flux [TO). A magnetic flux through the ring acts on Hilbert space states as \n){n'\ — > e'^"|n) (n'|e~'^" , and the 
corresponding "covariant derivative" governs the action (38 1. However, the Aharonov-Bohm analogy must not be 



carried too far: the operator [6, ] does not represent a genuine magnetic flux. Rather, it carries the information on 



the infinite extension of the "real" system, which is different from the single unit cell system for which (381 is defined. 
The difference to a real flux is that we are integrating over all gauge sectors (corresponding to field configurations 
Bq 0i for all values and 0') at once, and this restores the time reversal invariance of the system (which in the presence 
of a real flux would be lost.) These differences notwithstanding, the Aharonov-Bohm analogy will be a useful guiding 
principle in the discussion of the system below. 

Let us conclude this section with a general observation that will simplify our discussion below: in principle, the 
degrees of freedom of the theory are "infinite dimensional matrices", {Bg gi} depending on the continuous index 0. 



However in the absence of sources such as those specified in ( 30 ) , the integration over this matrix field gives unity. 



1. This is a consequence of supersymmetry. In the presence of a symmetry breaking source (the projector 



V in (30)), supersymmetry gets partially lost. Specifically, the decoupling of the source from all (^-indices other 



than 0± in (30 1 implies that the integration over the variables Bg_gi, {0,0') ^ {0+,0-) docs not contribute to the 



functional integral. These variables may simply be removed from the integration manifold (as long as we are probing a 



functional expectation value with fixed Bloch index configuration as in (30).) This reduction is a rigorous consequence 
of supersymmetry. In effect, it implies that the continuous index set {(0,0')} gets reduced to one element (0+,6'_), 
i.e. -Be, 61' — >■ Be^^g_Sg^e+6ei ^g_. Since 0± are fixed, we need not keep them as indices, i.e. Bne^a.ne^a' — > Bna.na'- 
Similarly, Qa0,a'9'{n) — >■ Qa,a'{n), where the reduced Q is generated by the reduced B. In the effective action (38), 
the reduced form of the operator reads 

= {0+E]^K + (^-Ef^) ® < (39) 

where the matrix represents the sign change in the Bloch "flux" under time reversal. Throughout, we will mostly 
work in the reduced representation. 



OFF RESONANCE 



In the next two sections, we will employ the effective field theory (38) to discuss the long-ranged physics of the 
system. Our discussion is organized in two parts: we first discuss the case of (near) irrational h/Air where ^/g 
and ^ Dq will be the localization length of the system. In this regime, periodicity effects relating to the q- 
commensurability of Planck's constant play no role and the system can be regarded as infinitely extended in angular 
momentum space. In a manner to be reviewed below, it then becomes similar to an infinitely extended disordered 
wire. This limit will serve as a benchmark against which to compare the physics of the g-periodic system. 

The commutators [0, ] appearing in the effective action have an eigenvalue spectrum limited by 2n/q. For ^/q 0, 
this is negligible in comparison to the excitation gap ~ of the localized systems, and we may simply ignore the 
^-contribution to the effective action. 



In this regime the effective action ( 38 ) reduces to the nonlinear cr-model action of quasi one-dimensional disordered 
metals (TS] [26]. To discuss the localization behavior of this system, we consider the density-density correlation 
function defined in ([s]). We first note that the system is "translationally invariant" in angular momentum space, 
meaning that Ki^{ni,n2) — K^^^ni — 112) depends only on coordinate differences. As such, it can be expressed in 



a Fourier representation, K^{(f)) = ^ '^'^"Ki^{n). Comparison with Eq. (10) shows that this bmction yields the 
energy increase as 
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By unitarity, the function K has to obey the Umit behavior hm0_>.o K^{(f>) 
compatible with this requirement (and time reversal symmetry) reads 



1 



—iuj + D{uj) 



The most general low (f> asymptotic 

(41) 



where D{ijj) is the so-called dynamical diffusion constant. 

Referring for a much more extensive discussion to Ref. |15| . we here briefly discuss the behavior of D{uj) upon 
lowering the frequency: 



• At high frequencies, ujDg 



oo, fluctuations of the fields Q are strongly damped and the effective action (38) 
may be approximated by the quadratic expansion around Q = cr^j^: 



S[B] 



dnsii{B{-Dqdl - iuj)B). 



Likewise, the correlation function ( 30 1 reduces to its Gaussian representation 

Ku:{ni,n2) = {stiiBin2)V)stiiB{ni)V)). 



(42) 



(43) 



Doing the Gaussian integral, we obtain (41 1 with D{lj) ~ Dq given by its bare value. Substitution into (40) 
shows that the kinetic energy increases as 



E{i) 



Dqt, 



(44) 



i.e. the variable E ^ tt? shows the behavior characteristic of diffusive dynamics. 



• Upon lowering w, but still keeping lo > ^, corrections to the bare value begin to form. The ensuing perturba- 
tion series is controlled by the parameter {ujDq)~^. Technically, these corrections may be obtained by expansion 
of the action ( 38 ) to higher order in the constituent fields B and subsequent integration against the Gaussian 



kernel (42). To leading order, this leads to a renormalization of the diffusion constant, 

1 



D{UJ) = Dq 



1 



27r —iu ■ 



(45) 



Physically, the correction term represents the modification of the diffusion constant by quantum weak localization 
corrections OIH]- The reduced diffusivity leads to a lowering of the energy increase, E{i) = Dqi— -j^^DqP^^ 

where i <^ Dg. 



Finally, at low frequencies, w ^ Dq ^, the function K^{(j)) has the form 

K^{<p) > : , 

Ai<l>) = 1 - C{3)Dy + 0{Dl<j^^), 
where C is the Riemann ^ function. It indicates that the diffusion constant approaches the value 



(46) 



(47) 



Eq. (47) implies vanishing diffusivity at low frequencies, a hallmark of strong dynamical localization. Substitu- 
tion into ( 40 ) indeed shows that 



Eii) D^q, 



(48) 



saturates at a constant value. At the saturation time ^ Dq, the two results (44) and (48) match. Therefore, 
the localization length scales as ^ ^ \J Dqij^ — Dq. 

To make the phenomenon of localization more explicit, one may consider the asymptotic form of the real time 
representation of the density correlator, K{n) = K{n, i — > oo). This function is given by [T51 
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where we have neglected the numerical prefactor and have defined 

e = DJ2 (50) 

as the localization length. This formula was first obtained by Shepelyansky [20j based on phenomenological 
arguments. 

Summarizing, the QKR supports diffusive excitations in angular momentum space at length scales An > K/h 3> 1 
before giving way to localization at scales An > ^ ^ Dq. These results hold for irrational values of h, and sufficiently 
large values of the kicking strength, K. For small values of Planck's constant, h <^ 1, the parameter Kg in (37) 

approaches the classical limit, Kg K, and the diffusion constant becomes 

^'^'(S) i^-^MK)^2jUK)+2Ji{K)), (51) 

where the corrections are the classical Rechester- White corrections [571 [2H]- In this limit, our result ^ = Dq/2 confirms 
an empirical formula by Shepelyansky and co-workers [2 [55V[5T| . 

VI. RESONANCES 

We now consider the system at rational h = Airp/q J2j. Pioneering work on the physics of the kicked rotor at these 
"resonant" values of Planck's constant has been done by IS [Hj (see also Ref. [5] for review.) To assist the orientation 
of the reader, we begin by reviewing some key results of that work in the language of the present formalism. We 
will then proceed to extend the theory so as to include the effects of diffusion and Anderson localization. (While the 
ramifications of localization at resonant values of h have been addressed before [SJ [7] , these discussions were based 
on phenomenology and numerical diagonalization.) Throughout this section, it will be convenient to represent the 



function E{t) as (cf. Eq. ([TO])) 

9-1 

2tt 

where K^{n) is the density correlation function evaluated at fixed parameters 0, 9' . 



E{i) = -l * ^ / {de)dll,^^KUn), (52) 



A. Theory of Izrailev and Shepelyansky 

In Ref. [6 , the periodicity effects in angular momentum space were discussed in terms of an eigenbasis of the 



effective Floquet operator (26). To reproduce these results in the language of the present formalism, we start from 
the formal diagonalization 

j 

where R{9) = {Rij{9)} are unitary q x q matrices parametrically depending on the Bloch angle, and £j{9) are the 



^-dependent eigenenergies. We may use the transformation R{9) to transform the fields Z in (31) to an "eigenbasis 
representation" , 

Zne,n'e' — >■ Rnj{(^)Zj9,j'e'R]fn'(^')- 

This transformation is unitary and does not affect the integration measure. As a result the Gaussian expansion of 
the action assumes the form 
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where we have included the frequency contribution to the quadratic action. Likewise, the Gaussian approximation to 



the density correlation function ( 30 1 is given by 



\de){de')e' 



where in the last step we did the Gaussian integral over the action S[Z]. This result is but a formal expansion of the 
density-density correlation function in a Bloch eigenbasis pT. It may be employed to obtain estimates for the energy 



increase E{t), as given by Eq. (10): 

This formula has been derived by IS [5]. Except for a few small values of q{~ 1,2,4), it cannot be evaluated in 
closed terms. However, the general structure of the result indicates that the energy stored in the system increases as 
in Eq. Q. Specifically, for large times the energy increases quadratically. To understand this asymptotic behavior, 
notice that at large times i 3> A~^, where the level spacing A ~ \ej — £j+i\ is a measure of the typical separation 



between neighboring levels, the dominant contribution to (53) comes from the diagonal term, j = j' . Describing the 
auto-correlation of individual levels (or "bands"), tj{0), this term is of "deep quantum origin". In this limit, the 
^-derivatives dominantly act on the rapidly oscillatory phases '-^ exp(iej(0')t), and we arrive at the estimate 

Notice that the sum X]n=o /fc('^^)l^"j'(^)l'^ \Ro,j{(^)\'^ is (approximately) normalized to unity. This means that 

^^i((9,(e,(^)))^), (54) 

can be interpreted as measure of the typical "level velocity" corresponding to the Bloch-phase dependence, and the 
average is taken over different quasi-energy levels (for fixed Bloch angle 9). Corrections to the dominantly quadratic 
increase arise from other derivative combinations, where the time derivative acts once (the linear term in (|4|) or twice 
(the term b{t)) on the diagonalizing matrices. 

However, at shorter times t < A^^, inter-band contributions j ^ j' begin to matter, and these can not be effectively 
described by a formal eigenfunction decomposition. In the next section, we present a n ana lytical theory of resonances 



in the regime i < q < ^ which yields the full time dependent profile E{t). In section VIC we will briefly address the 
regime of vary large periodicity intervals, q > S,. 

B. Resonances in the regime £ < q < ^ 



We will now address the physics of resonances within the framework of the theory ( 38 ) . Specifically, we will aim to 



obtain concrete expressions for the energy increase E(t) in terms of universal system parameters. 

1. Frequencies a; 2> A 

We start by considering the regime of frequencies a; ^ A. Recall that A = 27r/(7 is the effective quasi-energy level 
spacing (or the Bloch "band gap") of a periodicity volume. Thanks to the large value of uj, we may expand Q to 



quadratic order in B and consider the Gaussian approximation to the action, (42 1, where, however, the derivatives 



dn ^ dji + i[6, ] have to be generalized to covariant derivatives and the commutator acts as 
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Switching to a Fourier representation B{n) = q ^ J^tfi ^ and doing the Gaussian integral against the source 

K^{n)=q-^Y. [idO){d0') 



terms (43 1, we obtain for K^{n) = K^{n,Q) 



Dq{(l) + d-e'f 



We next substitute this result into Eq. (10) and perform the straightforward integrals over n and uj. As a result, we 



obtain Eq. ( 44 1 ; the periodicity of the system does not play a significant role here. 



2. Frequencies a; < A 

For w ^ A and flux parameters 0+ ~ 9- relevant to the computation of energy increase, the action of quadratic 
fluctuations is no longer small in comparison to unity. This signals a breakdown of perturbation theory. Indeed, we 



anticipate from our discussion of Section VIA that this is the regime, where the band-dispersion, dgej{6), of individual 
levels begins to matter. It is known jl5) that the fine structure of individual levels of hyperbolic systems is encoded 
in non-perturbative (in lu) fluctuations of the field theory. 

Technically, our setup bears similarity to the problem of "parametric correlations" [32]. There, one usually considers 
correlations in the density of states of a system at slightly different values of some control parameter. In the present 
context, the role of that parameter is taken by the Bloch phases 9±. (For earlier studies of Bloch-phase parameter 
correlations applied to the physics of periodically extended Hamiltonian chaotic systems, see Refs. [551 IM] .) 

To set the stage for the non-perturbative calculation, we first note that for w ^ A, inhomogeneous fluctuations 
(modes of finite momentum (p) have a typical action -Ec/A > 1 (the action of the first non-vanishing Perron- Frobenius 
eigenmode of a diffusive system). This is much larger than the action of the homogeneous mode, Q{n) = Q = const.. 



which means that inhomogeneous fluctuations can be neglected. Substitution of the constant configuration into (38) 
shows that the "zero-mode action" is given by 

S[Q] = ^striDg[9,Qf - 2iloQcjI^). 

Referring to Ref. [15] for an in-depth discussion of the technical details, we here merely note, that the integration 
over Q is best performed in a polar coordinate representation of the Q-matrices. The idea is to represent Q as 
Q = TZQi)R~^, where [7?., cr^p] = and the "radial" matrix Qo is given by 



cos 6 i sin Q 
— isinO cos 8 



(55) 



e 



AR 

711 \ a f ^ ^ \ a • / (71 (72 



6*22 / bp V ^ / t \02 Oi 



with < 6* < TT, 6*1^2 > 0. The action S[Q] is "rotationally invariant" in that it couples only to the radial coordinates 
(A = cos6', Ai_2 = cosh6'i_2) of the polar decomposition. Of these three variables, A e [—1,1] is "compact", while 
''^1,2 G [IjOo] parameterize the non-compact domain of the super integration manifold. 
Inserting the polar decomposition of Q into the zero-mode action, we find 

S = ^{D,[{M)\\l - A2) + {9+ + 9_f{\l - 1)] + 2iu:{\ - A1A2)}, 

where A6' = 6*+ — 9-. The action S then suggests that, for generic values of the phases 9± = Oi^Tijq) = ©(A), 
Q-fluctuations off-diagonal in T-space (characterized by A2) are penalized by a large action 0{DqA). Physically, 
this means that for fixed values of the phases 9± time reversal invariance is effectively broken (the symmetry gets 
restored only upon integration over all phases), and Cooperon-like fiuctuations (fluctuations off-diagonal in T-space) 
are strongly damped. Under these circumstances, we may restrict ourselves to T-diagonal fluctuations ("diffusons") 
and set A2 = 1. The action then collapses to 

s-^-^^i>^^-^')+^^ii^^^^h (56) 

and the density correlation function Ki^{A9) is reduced to 

KUA9) = ? r dAi dA^e-"^(^?-^^)+^^(^-^). (57) 
4 Ji J_i Ai — A 
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To proceed further we substitute this result into Eq. (52) and perform the remaining integral exactly. As a result, we 
obtain 



E{i) = 



2 ./ 27r a;2 

where D{uj) is given by (see Fig. [2j the top and middle panel) 

2 



1 



1 



(58) 



(59) 



This result for the frequency dependent diffusion constant in periodic structures was first obtained in Ref. |33| . Doing 
the integral over frequencies in Eq. (581 we obtain Eq. ([5| (see Fig. [2] bottom panel). 

Eq. ([5|^ is the main result of the present section. It accurately describes the function i?(i), up to corrections in 
exp(— i^ci) (the effect of fluctuations inhomogeneous across the unit cell). Eq. ^\ coincides (including the coefficients) 
with a conjecture on the mean square displacement of periodic quantum maps(mimicking one-dimensional quantum 
random walks) formulated by Wojcik and Dorfman |35j . That conjecture was based on the phenomenological modelling 
of quantum maps in terms of random matrix theory. Note that in the short time regime, x < 1 or i < g, the quantum 
correction to the leading (diffusive) terms scales as (tjqf' . As pointed out in Ref. |35j, this scaling reflects universal 
fluctuations of the quasi-energy spectra and eigenfunctions of the system. To the best of our knowledge, this effect 
has not been discussed in connection with the QKR. 



Turning to the long time regime, a; > 1 or t > (7, comparison with the formal result (54) leads to the identification 



2D„ 



(60) 



of the typical level velocity square. To understand the meaning of this equation, notice that the so-called Thouless 
conductance |36j of the periodicity volume 

can be understood in terms of the typical sensitivity of energy levels to changes in the boundary conditions V'('3') = 
'0(O)e*'*. In the present context (j> = q9 ^ d/A is set by the Bloch phase. Substituting this identification into the 
formula above and using that ~ ZJ^A^, we obtain (60). Qualitatively, our results may be extrapolated into the 
regime where £ is comparable to q. Even quantitatively, (60) is in good agreement with the numerical finding of 
Ref. [5], where {{deej{e)f) was found to be O.m"^ /{qh^). 



C. Resonances in the regime q > £, 



Let us now briefly explore what happens in regimes where the periodicity interval exceeds the localization length. 
Nominally, we are still dealing with a g-periodic quantum system whose spectrum can be organized in Bloch bands. 
However, the fact that wave functions overlap only exponentially weakly 0(exp(— g/(8^))) across the g- volume means 
that these bands are exponentially narrow [8j. Accordingly, the level velocities deej ^ exp(— 5/(8^)) are exponentially 
small meaning that the 0{P) increase of the energy in large volumes q ^ £, will be weighted by an exponentially 
small coefficient. We do not aim to analyze this correction any further. 

Turning to the dominant contributions to E{t), we expect a profile as in the g — > 00 system, diffusive spreading at 
times i < Dq followed by saturation at larger times. To see how this comes about, notice that the principal effect of 
the finiteness of g is a substitution </>—>■(/) + AO in the argument of the density correlation function K^^ (0) . As far 
as the function E{t) is concerned, this change plays no role, the reason being that the derivative in (40) needs to be 
performed at both (f) — and AO = 0. This shows that the function E{t) behaves as in the g — > 00 system, apart from 
the exponentially small tail corrections mentioned above. 



VII. CONCLUDING REMARKS 



Extending earlier work '14', we have introduced a microscopic and parametrically controlled theory of quantum 
interference and localization in the QKR. Emphasis has been put on the discussion of resonance phenomena occurring 
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FIG. 2: Real (top) and imaginary (middle) part of the dynamical diffusion constant at resonant values h — Anp/q. Bottom; 
crossover from linear to quadratic energy increase. The parameters are p = l,q — 100 (solid, in red) and p = 2,q — 201 
(dashed, in blue). In both cases, K = 5. 



at rational values h = Airp/q of Planck's constant. The prototypical theory [T^ did not distinguish between rational 
and irrational values of h/iir, a lack of resolution that has prompted criticism [9l[T6], and led to the suspicion that the 
theory was based on hidden ad hoc assumptions. The present generalization is formulated for rational configurations 
from the outset (the irrational limit can be realized in terms of large coprime values of q and p) . It is shovifn that the 
system bears similarity to an effectively disordered ring of circumference q and subject to an Aharonov-Bohm flux. 
(The integration over all values of the flux restores time reversal invariance.) We have obtained quantitative results 
for the rotor's kinetic energy E(t), an observable that carries detailed information on its time-dependent transport 
characteristics. For periodicities q smaller than the localization length ^ (but larger than the "transport mean free 
path" K/h), we have applied non-perturbative methods to express the full time dependent profile E[t) in terms of 
the universal scaling form Eq. ([5|. (Earlier results on Eit) ^ were formulated in terms of expansions of the long time 
asymptotics with largely unknown coefficients.) For large values q > S^, we have obtained results for the scaling of 
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E{t) that are quantitative up to corrections in exp(— g/^). The resulting profile conforms with general expectations 
on energy diffusion in a strongly localized system. 

In general, the theory above respects the relevant symmetries of the rotor (time reversal invariance, g-translational 
invariance, and the duality to near classical dynamics forming close to fundamental resonances at h = Airp). Correc- 
tions to the diffusion constant are resolved in agreement with earlier work [501 IMl 130] ■ At the lowest energy scales, 
the theory resolves the parametric dependence of individual levels in quantitative agreement with exact quantum 
mechanical calculations. In general, we believe that the criticism of the field theory approach to the rotor [3 [H] has 
been ill founded: the theory resolves microscopic structures at a resolution no "ad hoc" approach would achieve. 

It may be interesting to generalize the approach to other driven quantum systems, notably the two-sided kicked 
rotor ^Tf], the kicked Haper model f55], the double-kicked particle |3H], and the kicked top [ID] . 
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Appendix A: Derivation of Eq. ( |28[ ) 



For [M/,^^^] = we have str(VF™) = for odd m. Therefore, 2strln(l ±W) ^ strln(l - W'^). Applying this 
identity to S[Z]\^=o (cf. Eq. ([22])) we obtain 



S[Z]l=o 



--str< In 
2 I 



Z 
Z 



In 



U 



W \ z 



z 



(Al) 



Noticing 



we transform Eq. ( Al I to 



(gain + 1) = 



1 - 



Z 
Z 



C/t-1 



u-i 

' U+l 



Upon insertion of [/ = exp(— iG) this gets us to Eq. (28). 



Appendix B: Derivation of Eq. (30| 



We first note that the prototype correlation function ( 15 ) can be represented as 
where the source matrix X is defined by 



AR 



aE 



AR ' 



where, V^^ „2 = ^ni,n^n2,n is a projector onto angular momentum site n and E'^^ projects on the ft sector of the 
theory. Tracing the fate of this source contribution, we see that the action of the functional ( 19 1 picks up a contribution 



^^{X (g) E\^ +X^ (E) E^^)'$-^. 
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Further, the first of the logarithms in ( 20 ) gets modified according to 

ln(l - Zt-'^Z'^t) 
-> In (1 - (Z + E]^^ ® (a£;ii ® + aE^ ® V"')) 

x{Z + E^\, «) {aE]} ® 75"^ + aE^^ (g) V"^))) , 

where the identification Z = t^^Z^t has been used. At this point, we may carry out the derivatives with respect to 
a and a to obtain 

KUni,n2) =(str((l - ZZ)-^Z{Ei], (g) E^^ (E)V''')str{{l - ZZ)-^Z{E^\ (g) E^^ (E)V"')), 



where the angular brackets stand for functional averaging over the functional ( 22 ) , and we have been using the 
symmetry (21 1. In the invariant language of Q-matrices, this assumes the form 



1 



(str(g Ef^ (E) E^], (E) E].^ (g) 7^"^ )str(g E}^^ (g E^], g) E\^ g) V'')). 



Finally, we may resolve the projector matrices in the Bloch basis. Using Eq. (23 1, and noting that the Bloch- 
momentum is 'conserved', i.e. that the phases 0-entering the correlation function are pairwise equal, we thus obtain 
Eq. (pOl). 



Appendix C: Massive mode fluctuations 

In Section [TV] we have considered the action of "soft" fluctuations diagonal and long-ranged in angular momentum 
space. We will here extend the analysis of the action to more general fluctuations. The purpose of this discussion 
is twofold: First, we will show that massive mode fluctuations feed back into the action of soft fluctuations in that 
they change their diffusivity. Second, we need to explore the role of massive fluctuations by way of an a posteriori 
justiflcation of our identification of soft modes above. 

As we are going to show below, massive fluctuations are limited in angular momentum space to scales An ^ £ = 
1). This is smaller than the resonance periodicity volumes Q > ^ considered in this paper. To keep the notation 
simple, we will therefore neglect the Bloch structure of the fields Z throughout this appendix. In practice, this means 



that we will work with the basic representation Z = {Zatn,a't'n'} defined in (19) 



1. Massive modes 

We begin by introducing a few algebraic structures that will facilitate the analysis below. The fields Z = {Znn'} 
are matrices in the Hilbert space % of angular momentum states \n). As such, they may be viewed as elements of the 
Hilbert space A4 = H ^H* , where H* is the dual angular momentum space, i.e. the space of states The space 

A4 is spanned by the states, 

\n,n') = \n) (g {n'\. 
We will also need the dual of Ai, i.e. the space A^* of states 

(n, n'l = {n\ ® \n') . 

The natural scalar product and resolution of unity in M are given by, respectively 

{ni,n[\n2,n'2) = Sm^n2 Sn[,n'^_, 
^ |n,n')(n,n'| = 1^. (CI) 

nn' 
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Our analysis of fluctuations below will be conveniently performed in the Fourier transformed basis, 



in,n'\^ f (C2) 



Eq. (CI I implies the orthogonality and completeness relations 

(0, (t|</)', a') = 2^<5(0 
^ r ^\cb,a)i^,a\^lM- (C3) 

As discussed in the main text, modes of lowest action have zero angular momentum difference, (7 = 0, and wave 
numbers (j) < where (j)o is an upper cutoff to be self-consistently defined. To make the separation of "massive" and 
"massless" modes explicit, we introduce a projection operator ttq 



TTo^y ^|^,O)(0,O|, (C4) 

where d(f> = Jj^i^^^ d(l). It is straightforward to verify that this is a projection, i.e. ttq = tto. The projection on the 
complementary sector of "fast modes" is given hy n = 1^ — ttq, with the explicit representation 

^ = / ^l</',O)(0,O| + ^ / (C5) 

where d(j) = J^g<:\^\<:.^ and J dip = J^^ d(f>. The projection operation defines a decomposition 

Z = TToZ + 7tZ = B + C, (C6) 
where B = ttqZ and C = ttZ are the massless and massive contribution to Z, respectively. 

2. Renormalized Gaussian action 

With these structures in store, we will now re-iterate the derivation of the Gaussian action. However, our discussion 
will extend that of Section IV A in that we will keep explicit track of C-fluctuations. Ignoring contributions 0{ujd'^) 
(which for the frequencies cj ^ 1 of physical interest are of no relevance), the quadratic action is given by S[Z] « 
Sfi[Z] + S^[Z] with 

Sfi[Z] = ^stT{ZZ - ZU''ZU) , (C7) 
^^[Z] = -yStr(ZZ). (C8) 

The action S^j — We first notice that the projector property of tt implies a decomposition 

S^[Z] ^ S^[B] + S^[C]. 
(The proof of this relation is based on the orthogonality property 

str(Z7rZ) = str{ZTr^Z) = str(7rZ7rZ), (C9) 
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which in turn fohows from the definition (C4).) In principle, the fast field contribution S'cj[C] leads to a modification of 
the fast fluctuation action to be discussed below. For the small values of w ^ Dq^ relevant to our present discussion, 
however, this modification is irrelevant. We therefore dispose of the contribution S^[C] and turn to the slow action 
Su: [B] . Expanding the fields B as 



B = 



d<j) 
2^ 



S(<^)|0,O), 



(CIO) 



where B{(j)) = (0, Q\B — {(f>, 0\Z and using the auxiliary identity (a direct consequence of Eq. (CI I and the definition 



tr(|0,a)Z) = (-0,-a|Z)e-^*° 



we obtain 



S^[B] 



iuj f d4> 



str{Bi^)Bi-(b)). 



2 J, 2TT 

At this point, we may transform back to angular momentum space. Abbreviating 

B{n) = Bn,n, 

and approximating momentum sums by integrations, ^„ — )■ / dn, this gives 



SJB]^-- 



dn stiiB {n)B{n)). 



(ClI) 



The action Sa — We next turn to the more involved discussion of the fluctuation action. Introducing the notation 



Ad^X = O^XO, the latter may be written 



as 



Sn[B, C] = -str [{B + C) {l - Ad^) {B + C) 

= lstr[B (l - Adj^) B - CnAdjjB - 7rAd^_iBC + C* (l - TrAd^) C] 



(C12) 



where the second line is based on ( C9 1 . 



We next integrate out the massive modes to obtain an effective action Sa [B] , 



(C13) 



The Gaussian integration over C can be performed exactly. To fulfill this task, it is convenient to use the saddle point 
method. The solution of the saddle point equation 



SSfi[B,C] 
5C 



-nMu-iB+{l - TrAdjy-OC = 0, 



is given by 



Cci = (I - irAdfjA 7rAdr7_iB. 



Applying Eq. (21) to this equation, we obtain 

Cci = (l-7rAd^)"VAdj;B. 

As usual in supersymmetric theories the integration over quadratic fluctuations around the saddle point gives a 
factor of unity, i.e. the eff ective action is obtained by straightforward substitution of the saddle point configurations 
(Z™, Z^) into the action (C12). As a result, we obtain 



1 



S'fl[B] = -str \B\l-kdfj 



1 



2"" V V "^^l-7rAd/ 



B 
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In the following, it will be convenient to split the action into two contributions, 



Sn[B]=S°a[B] + Sm: (C14) 
Sl!'\B] - -istr (b( Ad,y^^^^ 1 B 



where the two terms represent the native slow mode action {S^) and its renormalization by the inclusion of massive 



modes (S'^'), respectively. The structure of the action (C14| solely relies on the projector property of tt, but not on 
details of the projection operation. We next aim to bring the action into a more explicit form, and to this end, we 
need to take the dynamical structure of the Floquet operator into account. 

Before turning to the calculation, it may be worthwhile to discuss the physical meaning of the action 5^". This 
action describes the renormalization of the soft mode action by "massive" fluctuations. Consider a generic massive 



fluctuation mode Cmn^, non-diagonal in angular momentum space. As discussed in Section IV A we may loosely 
interpret Cnin2 the representative of quantum states |ni)(n2|. Passing to the Wigner representation, the off- 
diagonality, i.e.. An = ni — n2 7^ 0, leads to an excitation in the angular direction of the "phase space". As we 
will see, the kernel (1 — TrAd^y)"^ describes the ways in which the memory of initial "phase space" information 
is maintained. For generic irrational values h <C these effects diminish rapidly. (One may construct certain 
infinite fraction representations of h which sum up to irrational values and, nonetheless, lead to long-ranged angular 
correlations |23| . At this point, we are not able to get the ensuing memory kernels under control.) The inclusion 
of this memory kernel will lead to a renormalization of the diffusion constant of the theory. Equally important, it 
may affect the length scales at which the dynamics becomes genuinely diffusive. Again, the situation bears analogy 
to the physics of disordered metals. There, wave packets initially centered around a definite momentum maintain 
the memory of their direction of propagation over a certain time scale, the so-called transport scattering time. In a 
manner largely analogous to our present discussion, the transport scattering time can be seen to determine an effective 
diffusion constant. 

The action Sf^[B] — We first consider the unrenormalized massless action 5". The derivation of this action 
essentially repeats the construction of Section |IV A| in the present notation and introduces technical building blocks 
relevant to the discussion below. Specifically, we will need the matrix elements of the adjoint Floquet operator in 
the two bases {\n,n')} and {|0,(t)}. Using that {n+,n-\kA^\n'_^,n'_) = {nj^\lJ\n'j^){n'_\ij^n-) , a straightforward if 
somewhat tedious calculation obtains 

'K\ ( K 

^n'_ —n 



(</.,a|Ad^|0',a') = 2^^(0-0'-f /ia)e-*§'^'+^T(--')j^_^, (^^sin(^^^^ . (C15) 



where 5(x) stands for that the (5-function holds only modulo 2tt, and 



are the Bessel functions of the first kind. Using this result and representing the slow fields as (CIO I, we then obtain 



Sl[B] = \ I ^^str(i3(0W))(-0,O|l-Ad^|0',O). 



Using Eqs. (C3) and (C15l this becomes 



-ll^MBim^cb)) (1- Jo (^sin(^-0)) . (CI7) 

This action will describe diffusive correlations of the field B, provided the angular kernel (1 — Jq) can be approximated 
as a quadratic function of cf). This is the case, if </> is sufficiently small. Specifically, for a soft mode angular cutoff 

'^o«ri = -|, (CI8) 
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with h/K sufRciently small, we have ^| sin ^| ^ 1 which means that the Bessel function may be expanded as 



to arrive at Sf^[Z] ~ J ^0^str(i3(0)i3(— 0)). Finally, transforming back to angular momentum space, we obtain 



Dn 



dn siY{dnBdnB) 



(C19) 



Notice that the condition (C18l stabilizing this two-derivative action has an analog in the physics of disordered 



metals. There, the wave numbers of diffusive fluctuations (fluctuations governed by a second order derivative evolution 
operator) are smaller than the inverse mean free path, /f^? due to impurity scattering. In the present context, the 
scale K/h plays a role analogous to Itr and modes fluctuating at larger length scales evolve diffusively. 



The action Sf^^[B] — We next discuss the second contribution Sl^ to the fluctuation action (C14). Inserting the 



expansion B = 
starting point 



/s '^■^i4>)\'t'^^) into Eq. (C14| and expanding the resolvent operator (1 



^Ad^) 



we obtain our 



5" = -: 



27r 27r 



str(B(-<^)B(0')) £(</), 0|Ad^(7rAdj>)"|<^',0). 



n=l 



We now insert resolutions of unity ( C3 1 between the n factors TrAdj^, use that the projectors tt appearing are diagonal 
in the basis a)}, i.e. 7r|0, cr) = 7r(0, cr) |0, cr), where 7r((/), cr) — (5ct,o6'(|</>| — </'o) + (1 — 5^^), {0{x) is the Hea viside 
fimction, which is unity for a; > and zero otherwise.) and substitute the explicit form of the matrix elements (C15) 
to arrive at the "path integral" representation 



27r 27r 



str(i3(~0)S(^')) 



n+l 



2K . 

sm 

h 



,(C20) 



Here, the integration measure is given by 



and the boundary conditions (0„+i,cr„+i) ~ (0,0) and {(j)o,ao) = (0',O) are understood. We may now use the 
(5-function constraint in the integral to collapse the phase factor to multiple 27r, 

n+l 



haf 



Omod27r. 



The constraint further implies 



t>' — h CTj mod 2tt. 



(C21) 



For generic irrational /i/(47r) (for "non-generic" values, see Ref. [23]), is not close to a multiple 27r which implies 
^ sin(0.,/2) = 0{j) > 1. Now for large arguments. 



|a;|>l 1 



which means that Eq. ( C20 ) is a controlled expansion in {h/Kf/'^ < 1. We may truncate this expansion at any 



desired order n. In conjunction with the near irrationality of h/Air, the finiteness of the sum X]J=i '^j ^ N means that 
(jiJ2'j=i ^j) niod 27r, can not be of 0{4>') < 0(h/K) unless X]j=i — 0- We thus conclude 



(C22) 



25 



giving 4> = 4>''- the integral over intermediate states conserves the momentum <j). 
The considerations above fix the simplified action 



S^[Z]=-\ I ^str(i?(-0)B(0))f;^n J..-.._, (^sin 

n=l {a} i=l ^ 



(C23) 



where (pi is given by (C21) and J2{a} ~ Hi E<j j subject to the constraint (C22) 



1- 


• 1- 


^1, 




-1- 


i \^2^^/^ i 


-2- 





FIG. 3: Path of lowest order contributing to the renormalization of the diffusion constant. 
It is instructive, to think of the set {{4>i, Ci)} in terms of a path 
((/),0) ^ ((/)- h<7i,<Ji) ^ (0- hai - ha2,(T2) > {(p - hai /ict„_i, cr„_i) ((/',cr„) (0,0). (C24) 

We may imagine the n + 1 links "— ?>" of this path to be weighted by Bessel functions. Specifically, the first and last 
link, respectively carry weights J^^ ^^sin^|^^ and J^a„ (x^^i'^(f))- These weights imply a further effective 
constraint on the sum over configurations {tr}: since (j) < 4>o <ti h/ K, the terminating Bessel functions carry arguments 



2K I 



sin 1 1 ~ \K(f>/h\ 1. We may thus apply the asymptotic formula 



|eKi e 



kl 



2l'^lcr! 

to conclude that the dominant contribution to the sum has 



The terminal weights then yield a factor 



2K . 
sm , 
h V 2 



kil = |ct„| = 1. 



2K . 

sm , 
h V 2 



(C25) 



V 2h 



sgn((Ti cr„). 



For the arguments of the inner Bessel functions, regardless of the value of h, they much exceed unity, which means 
that the sum above will be dominated by the shortest paths obeying the criteria ( C22[ ) and ( |C25 ). 

The contributions of leading and sub-leading order to the path sum are shown in Fig. [3j Evaluating these paths in 
(|C23|), we obtain 



sm = 



Dn 



d4 
2^ 



02 str(B(-0)B(0)) {-2J2{Kq) - 2Jl[Kq) + 2Jl{Kq)) 



(C26) 



Passing back to angular momentum space and adding to S^^ the compounds (C19l and (Cll), we arrive at the final 
form of the Gaussian action 



S[B] 



dnstT{B{n){~Dgdl ~ iuj)B{n)), 



where the diffusion constant is given by ( 36 ). Corrections to the above diffusion constant can be computed by including 
path of higher complexity and will be of 0{Kq ^^'^). However, at this order in the expansion, anharmonic fiuctuations 
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of 0{sti{BCBC)) begin to play a role and the computation becomes significantly more complicated (see Ref. [20] for 
the form of the ensuing corrections.) 
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